Optomechanical transduction of an integrated silicon cantilever probe using a microdisk resonator 
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Sensitive transduction of the motion of a microscale cantilever is central to many applications in mass, 
force, magnetic resonance, and displacement sensing. Reducing cantilever size to nanoscale dimensions 
can improve the bandwidth and sensitivity of techniques like atomic force microscopy, but current optical 
transduction methods suffer when the cantilever is small compared to the achievable spot size. Here, 
we demonstrate sensitive optical transduction in a monolithic cavity-optomechanical system in which a 
sub-picogram silicon cantilever with a sharp probe tip is separated from a microdisk optical resonator 
by a nanoscale gap. High quality factor (Q « 10 5 ) microdisk optical modes transduce the cantilever's 
MHz frequency thermally-driven vibrations with a displacement sensitivity of « 4.4 x 10 -16 m/\/Hz and 
bandwidth > 1 GHz, and a dynamic range > 10 6 is estimated for a 1 s measurement. Optically-induced 
stiffening due to the strong optomechanical interaction is observed, and engineering of probe dynamics 
through cantilever design and electrostatic actuation is illustrated. 



Micro- and nanoscale cantilevers are at the heart of many 
applications in mass, force, magnetic resonance, and displace- 
ment sensin g 1 ^ 4 . In atomic force microscopy (AFM) 5 , the 
push towards smaller cantilevers^ is motivated by the abil- 
ity to increase mechanical frequencies while maintaining a 
desired level of stiffness. This influences the force sensi- 
tivity and measurement bandwidth, in turn determining the 
image acquisition rate and ability to resolve time-dependent 
forces and acquire additional information about the tip- sample 
interaction potential 8 . Standard optical methods for trans- 
ducing cantilever motion include beam deflection 9 and laser 
interferometryiS, and in macroscopic devices that are 1 mm 
x 1 mm x 60 jum (length, width, and height), quantum- 
limited displacement sensitivity of 4x 10 -19 m/^/Hz has been 
achieved^. Interferometric approaches using a high numer- 
ical aperture objective have also been used in micro- scale 
devices, resulting in displacement sensitivities of 3xl0 -14 
m/\/Hz for cantilevers that are 20 jum x 4 ^mx 0.2 jum and 
1 x 10 -15 ui/a/Hz for larger conventional cantilevers (223 ^m 
x 31 ^mx 6.7 ^m) 12 . However, as the cantilever dimensions 
are pushed below the detection wavelength, diffraction effects 
limit the sensitivity of these approaches^, and near-field op- 
tics and/or integrated on-chip detection methods can be of sig- 
nificant benefit. 

To that end, researchers have recently used evanes- 
cently coupled on-chip waveguides 14 acting as doubly- 
clamped cantilever s 15 ! 16 to demonstrate displacement sensi- 
tivities of 3.5x10" 14 m/V^Hz, while end-to-end waveguides 
acting as singly-clamped devices 17 have achieved similar 
performance 18 . Although these waveguide-based approaches 
are optically broadband, the strong, multi-pass interaction 
provided by optical cavities can be of considerable advan- 
tage. Cavity optomechanic a 1 9 ^ 2 1 has seen substantial re- 
cent progress, where in many cases the optical resonator also 
acts as a mechanical oscillator, and its internal vibrations 
have been transduced with measurement imprecision at or be- 
low the standard quantum limit 22,23 and with absolute dis- 
placement sensitivities in the 10 -17 m/\/Hz to 10 -18 m/V^Hz 
range 24,25 . In contrast, here we focus on transducing the mo- 
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FIG. 1: (a) Scanning electron micrographs of the cantilever- 
microdisk system. The right image has the FEM-calculated z- 
component of the magnetic field for the TEi 51 mode overlaid on the 
structure, while the inset shows a zoomed-in region with cantilever 
width W and gap G; (b) Simulated mechanical modes (amplitude ex- 
aggerated for clarity) with dominant displacement along the x-axis 
for W=65 nm; (c) Predicted optomechanical coupling gQM between 
the h\ x cantilever mode and TE/TM modes of the microdisk. 



tion of a cantilever probe, requiring a design in which the 
cantilever can be brought near a surface and its fluctuations 
sensed by a nearby optical cavity without inducing excessive 
optical loss. 

A similar approach was presented in Ref. 26| where 
doubly-clamped SiN x nanobeams were brought into the near- 
field of SiC>2 microtoroid cavities fabricated on a separate 
chip. In comparison, here we fabricate a cantilever-optical 
cavity system on a single silicon device layer, while tai- 
loring the cantilever geometry for both strong optomechan- 
ical interactions and applicability to AFM. Beyond demon- 



2 



TABLE I: Calculated and measured properties of the h mx cantilever modes. 



Mode k (calc.) m (calc.) £2m/271 (calc.) (expt.) Tm/2tz (expt.) Qm (expt.) 





0.14 N/m 


0.73 pg 


2.23 MHz 


2.35 MHz 


479 ±8 kHz 


4.9 


hx 


1.41 N/m 


0.58 pg 


7.82 MHz 


7.89 MHz 


598±11 kHz 


13.1 


hy x 


5.72 N/m 


0.35 pg 


20.37 MHz 


20.51 MHz 


533 ±2 kHz 


38.5 


Kx 


12.43 N/m 


0.32 pg 


31.17 MHz 


31.36 MHz 


706 ±4 kHz 


44.4 


hsx 


41.95 N/m 


0.44 pg 


49.36 MHz 


49.86 MHz 


815±4kHz 


61.2 


hx 


74.05 N/m 


0.40 pg 


68.13 MHz 


68.71 MHz 


752 ±43 kHz 


91.0 



strating sub-fm/V^lz sensitivity to cantilever motion, this ap- 
proach has many potential benefits for AFM. Silicon's high 
refractive index allows for significantly smaller optical cav- 
ities to be used, yielding stronger cantilever-cavity coupling 
rates and permitting higher bandwidth operation. Moving 
to silicon opens up potentially advanced device functional- 
ity, including electrostatic actuation and integrated optical 
waveguide readout. By largely separating the mechanical 
and optical designs, engineering of the cantilever geometry 
to achieve desired parameters can be accomplished without 
adversely affecting the optical readout mechanism. In addi- 
tion, the strong optomechanical interaction can allow for op- 
tical control of cantilever mechanics, through effects such as 
optically-induced stiffening and optically-driven mechanical 
vibration s 19 i 2Q i 24 i 27 . Finally, this platform provides simplifi- 
cations with respect to free-space detection systems that may 
improve measurement stability and be of importance in par- 
allelized multi-probe measurements^ or environments with 
limited optical access. This work lays the foundations for a 
class of practical nanoscale mechanical sensors enabled by 
cavity optomechanics. 

Device geometry and simulation 

A simple device geometry is shown in Fig. [Ha), with fabri- 
cation details given in the Methods. A semicircular cantilever 
of width W is suspended at its ends and separated by a gap G 
from a 10 jum diameter silicon microdisk. The silicon is 260 
nm thick, and the cantilever has been designed to support a 
sharp tip at its midpoint. Devices are fabricated with W=65 
nm, 100 nm, and 200 nm, and nominal values G=50 nm, 75 
nm, and 100 nm. Scanning electron microscope (SEM) im- 
ages indicate that W is typically within ±5 nm of its nomi- 
nal value, while G is often smaller than the nominal value by 
a couple tens of nanometers, though charging effects due to 
the electron beam limit this estimate. The cantilever geome- 
try is chosen to maximize its interaction with microdisk opti- 
cal modes while minimizing the scattering loss induced by its 
presence. Optical modes are labeled TE PiH and TM p ^ n , accord- 
ing to polarization (transverse electric or transverse magnetic) 
and radial (p) and azimuthal (n) order. Three-dimensional fi- 
nite element method (FEM) eigenfrequency simulations indi- 
cate that, for W=65 nm or W=100 nm, cavity quality factors 
(2s) in excess of 10 6 can be achieved for TEi jW and TE2, n 
modes, and Qs in excess of 10 5 can be achieved for TMi 5 „ 
modes, even as G decreases to ^30 nm. In comparison and as 



a baseline, fabricated microdisks without cantilevers exhibit 
Qs in the mid-10 5 to low-10 6 range. 

Mechanical modes of a W=65 nm cantilever are deter- 
mined from FEM simulations (see Methods), with represen- 
tative modes shown in Fig. \V[b). We have focused on the h mx 
modes, which are even symmetry in-plane modes whose pri- 
mary displacement direction is normal to the gap (x direction), 
as they are the dominant modes that are optically transduced 
and are of particular relevance to AFM work. The predicted 
stiffness (k), resonant frequency (fljjf), and effective mass (m) 
of these modes are compiled in Table U Focusing on the h\ x 
mode at = 2.23 MHz, its optomechanical coupling to 

p=l and p=2 optical modes in the 1550 nm band, defined as 
gOM = d(ti c /dG (co c is the cavity mode frequency), is calcu- 
lated by FEM simulation and displayed in Fig. QIc). For the 
range of gaps studied in this work, £om/27C & 0.5 GHz/nm 
to £om/271 ~ 3.0 GHz/nm. This is about two orders of mag- 
nitude larger than #om for SiN x cantilevers coupled to SiC>2 
micro toroids 26 , and is due to the more tightly confined optical 
modes supported by the silicon microdisks. 
Transduction of cantilever motion 

We measure the fabricated devices using a fiber taper cou- 
pling method 2 - shown schematically in Fig. Oa) (see Meth- 
ods). A 1550 nm band tunable diode laser is attenuated and 
coupled into the devices using an optical fiber taper waveg- 
uide, a single mode optical fiber whose minimum diameter has 
been adiabatically and symmetrically reduced to around 1 jam. 
At this diameter, the waveguide mode's spatial profile extends 
well beyond the glass core into the surrounding air cladding, 
and this evanescent tail is used to excite and collect light from 
the microdisk modes. The signal exiting the cavity is split by 
a 90:10 fiber coupler, with 10 % of the light used for moni- 
toring the transmission level and recording swept- wavelength 
transmission spectra, and 90 % sent into a radio frequency 
(RF) photodetector, after which an electronic spectrum ana- 
lyzer measures RF oscillations in the detected signal. 

Normalized transmission spectra over the full wavelength 
band for TE and TM polarized modes of a W = 65 nm, G = 50 
nm device are shown in Fig. [2b), along with zoomed-in scans 
of individual modes. The polarization of the modes is de- 
termined by comparing the free spectral ranges for modes 
of a given radial order with those predicted from simula- 
tion. Loaded cavity Qs of 8.0 x 10 4 and 1.8 x 10 5 are observed 
for this device (corresponding intrinsic Qs of 1.1 xlO 5 and 
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FIG. 2: (a) Setup for device characterization, (b) Broad wavelength scan (left) for TE (top) and TM (bottom) modes of a typical disk-cantilever 
device (W=65 nm, G=50 nm). Zoomed-in scans (right) show data (green) along with a doublet model fit (black), (c) Broad RF spectrum of 
a disk-cantilever device (W=65 nm, G=100 nm), transduced by fixing the probe laser on the short wavelength side of the TE-polarized mode 
shown in the inset (black=low power, Pi n = 14. 1/^ W, red=high power, Pi n = 223/j W). Mechanical modes below 100 MHz (blue) are due to 
the cantilever, while modes at 364.63 MHz and 577.20 MHz (green) are due to the disk, (d) RF spectrum of a disk without the cantilever, 
displaying modes at 364.74 MHz and 576.20 MHz. The inset shows a high-<2 TE optical mode of the disk (blue) with fit (black), (e) Zoomed-in 
RF spectrum of the disk-cantilever, showing the h mx modes (blue), calibration peak (purple), and detection background (black). 



2.1 xlO 5 , respectively), which supports doublet modes due 
to surface-roughness-induced backscattering that couples the 
clockwise and counterclockwise modes of the cavity 30 . Over 
all devices, Qs of 5xl0 4 to 2xl0 5 are typically observed 
for TEi jW , TE2,rc and TM\ iH modes, though occasional devices 
have Qs as high as w 6x 10 5 (see supplemental data). Optical 
transduction of the cantilever's motion due to thermal noise is 
performed by fixing the laser on the blue-detuned shoulder of 
a TE-polarized cavity mode. A 1 MHz to 600 MHz spectrum 
for a W = 65 nm, G = 100 nm device is shown in Fig. [2fc), 
and contains several peaks. Those below 100 MHz originate 
from motion of the cantilever, while those at higher frequen- 
cies (364.63 ±0.35 MHz and 577.20 ±0.25 MHz) are from 
motion of the disk. This is confirmed by measuring the RF 
spectrum of a disk without a cantilever (Fig. [2d)) through a 
high-g cavity mode (loaded g = 5.7xl0 5 ±0.5xl0 5 , intrin- 
sic Q w 1 .0 x 10 6 ), which yields RF peaks at near-identical fre- 
quencies (364.74 ± 0.03 MHz and 576.20 ± 0.03 MHz). FEM 
simulations indicate that the higher frequency mode is the 
disk's radial breathing mode (RBM); its measured linewidth 
is Vm/Ik = 21.68 ± 0.06 MHz, corresponding to Qm ~ 27. 

Focusing on the frequency range between 100 kHz and 100 
MHz, a higher resolution RF spectrum at 223 juW of input 
power (Pin) into the cavity is shown in Fig. [2e). The frequen- 



cies of the transduced modes (Table UJ) correspond well with 
the previously described simulation results. The mechanical 
quality factors of these modes are between Qm ~ 5 for the h\ x 
mode and Qm ~ 61 for the hs x mode (Tabled]); these values are 
likely limited by air damping^! The detection background, 
shown in Fig. [2e) in black, is found by placing the laser 
off-resonance while maintaining a fixed detected power. Fo- 
cusing on the h\ x mode, its calculated effective mass and mea- 
sured freque ncy cor respond to a peak displacement amplitude 
of 

Xrms — ykfiT/k ~ 160 pm when driven by thermal noise 
at 300 K. We use x rms and Tm to convert the RF amplitude 
in Fig. [2e) to displacement sensitivity 1 . The corresponding 
photodetector-limited sensitivity is 4.4 x 10~ 16 ± 0.3 x 10~ 16 
m/V^Hz. This value is consistent with that determined by a 
phase modulator calibration (Methods) to within our uncer- 
tainty in the disk-cantilever gap. It represents an improvement 
by about a factor of 100 with respect to other on-chip silicon 
cantilever experiment al 18 , is at the same absolutely sensi- 
tivity level demonstrated for SiN x cantilevers transduced by 
silica microtoroids^, and is about a factor of 5 times larger 
than the standard quantum limit 3 - for our system. Along with 
the sensitivity, two other important quantities that character- 
ize this system for its use as a displacement sensor are its dy- 
namic range and bandwidth. The maximum detectable dis- 
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FIG. 3: (a) RF spectra from a device (W=65 nm, G=75 nm) with 
Pin — 446 jjW at different laser-cavity detunings AX. (b) Image plot 
of the RF spectra as a function of AX. The cavity mode used for 
transduction is shown in the inset, (c) Zoomed-in portion of the im- 
age plot for the h\ x mode, showing optically-induced stiffening, (d) 
Zoomed-in image plot for the h\ x mode of a W=65 nm, G=50 nm de- 
vice. In (c)-(d), the RF spectra are displayed on a linear scale, with 
each spectrum normalized to the peak amplitude for that value of AA,. 



placement is approximately the ratio of the cavity linewidth 
(T/2% = 2.44 GHz) to goM, and is « 4 nm, giving a dynamic 
range > 10 6 (60 dB) for a Is measurement. The bandwidth 
(BW) is limited by the cavity's response time, which deter- 
mines how quickly it can transduce mechanical motion. We 
therefore expect a BW> 1 GHz, and this is substantiated by 
transduction of the 575 MHz oscillations of the disk as pre- 
viously described in Fig. Oc)-(d). Adjusting the BW (e.g., 
through the waveguide coupling) allows for gain/BW trade- 
off within the fixed gain-BW product. The large BW of these 
devices is one advantage of relying on large #om rather than 
ultra-high- Q for displacement detection. 
Optically-induced stiffening 

Increasing the optical power coupled into the cavity causes 
several notable changes in the RF spectrum, as seen in Fig. 
Oa)-(b) for a W=65 nm, G = 75 nm device, where the cou- 
pled power is changed by fixing P[ n = 446 juW and varying 
the detuning AX between the laser and cavity mode. First, the 
spectral position of the h\ x mode changes from £Im/2k ^2.24 
MHz at large A A, to Qf M /2n « 3.26 MHz at AX = -21 pm be- 
fore returning to close to its original value at near-zero AX 
(Fig. [2c)). One explanation for this is the optical spring ef- 
fect, an optically-generated rigidity of the mechanical oscilla- 
tor, as seen in other work o 19 ^ 27 ! 33 . In particular, if we take 
the measured values for £Im, AA,, T, and internal cavity 
energy U (determined by P{ n , transmission contrast, and F), 
the value of #om that best matches the maximum frequency 
shift is goM /27T = 1 .4 GHz/nm, corresponding to a gap G « 60 
nm for the TE2 7 45 mode. Similarly, Fig. [3td) shows a shift 
from Q m /2k w 2.26 MHz at large AX to Ol M /2% « 4.25 MHz 



at AX = —41 pm, in this case for the h\ x mode of a W=65 nm, 
G =50 nm device. This shift is consistent with ^om/2ti = 3.0 
GHz/nm, corresponding to a gap G « 32 nm for the TE2 7 45 
mode. Both of these gaps are smaller than the nominal val- 
ues, but are reasonable given the variation observed in SEM 
images of fabricated devices. 

Along with the change in frequency, the linewidth of the 
h\ x mode changes from Tm/271 ~ 410 kHz at AX = —61 pm 
to Tm/2tz w 860 kHz at AX = —21 pm, indicating damping. In 
addition, the increase in RF amplitude of the h mx modes is ac- 
companied by a broad background which, in certain detuning 
ranges, produces peaks in the RF spectrum not seen at lower 
powers and at frequencies that are not predicted by mechani- 
cal simulations of the cantilever. The precise nature of these 
effects is not understood, though a likely cause is the inter- 
play between free-carrier and thermal effects that takes place 
in silicon microdisks as the intracavity energy is increased. 
Measurements of devices with and without cantilevers (sup- 
plementary information) show behavior consistent with previ- 
ous observation of such effects 34 . It should also be noted that 
thermal effects have been observed to generate damping for 
blue-detuned excitation in other optomechanical systems 24 . 
Cantilever engineering and outlook 

While optically-induced stiffening provides real-time con- 
trol of the cantilever properties over a certain range, a number 
of modifications to its geometry can improve its applicability 
to different AFM applications. The sub-N/m spring constant 
of the h\ x mode is suitable for weak force measurements in 
which the cantilever is undriven, but in dynamic techniques 
for which the best imaging conditions have been achieved, 
such as frequency modulation AFM^, spring constants in the 
tens of N/m to hundreds of N/m range (or more) are desir- 
able for small amplitude operation 5 . In our geometry, the can- 
tilever stiffness may be increased by increasing its width; figs. 
Eta)-(b) show the mechanical mode spectra for W=100 nm 
and W=200 nm devices. The h\ x modes at £Im/2% = 3.33 
MHz and £Im/2% = 6.96 MHz agree well with the simulated 
values of 3.42 MHz and 7.22 MHz. Based on the calculated 
effective masses, these values correspond to a cantilever stiff- 
ness of 0.52 N/m and 4.11 N/m, respectively, with the latter 
being a 30 x increase in stiffness relative to the h\ x mode of 
the W = 65 nm device. Stiffer cantilevers can be produced 
by a further increase in W, though degradation in the optical 
Q is expected unless G is increased, which can then limit the 
displacement sensitivity due to a reduced #om- Another op- 
tion is to use smaller diameter microdisks, to reduce the can- 
tilever length between its suspension points. Bare microdisks 
have radiation-limited Qs > 10 6 until their diameters are just a 
couple of micrometers 36 , and simulations predict that the h\ x 
mode of a W=100 nm cantilever coupled to a 4.5 jum diame- 
ter disk will occur at 7.96 MHz (k = 1.8 N/m). Another im- 
portant consideration is the modal structure of the cantilever. 
Though we have focused on the h\ x mode due to its displace- 
ment profile and transduction under thermal noise, in an AFM 
setting, the cantilever motion will be defined by both its ac- 
tuation mechanism and the surface it is interrogating, and its 
motion will be a superposition of its modes. This includes 
out-of-plane (z direction) and orthogonal in-plane (y direc- 



5 




5 10 15 

Frequency (MHz) 



FIG. 4: (a) RF spectrum from a disk-cantilever with W=100 nm, 
G=200 nm. Inset shows displacement profiles for cantilever modes 
not shown in Fig. \T\b). (b) RF spectrum from a device with W=200 
nm, G=200 nm. (c) RF spectrum from a disk-double-cantilever with 
W=65 nm, G=50 nm. Inset shows SEMs of the device geometry. 



tion) modes such as those seen weakly in the RF spectra of 
Fig. Hta)-(b). Engineering of the cantilever support geome- 
try can better isolate the h\ x mode in frequency space. The 
double cantilever structure shown in the SEM images of Fig. 
HJc) has h\ x as its lowest frequency mode, with the first out- 
of-plane mode vi significantly stiffened and shifted to higher 
frequencies. Figure HJc) shows an optically-transduced RF 
spectrum for such a device, with W=65 nm. Going forward, 
further modifications may be made to increase the stiffness of 
the cantilever, for example, through multiple short supports to 
surrounding areas. By combining this approach with W=200 
nm cantilevers and/or smaller diameter disks, we expect that 
&=100 N/m devices will be achievable. On the opposite end of 
the spectrum, very soft cantilevers are also of considerable in- 
terest, due to their application in measurements of very small 
forces 37 such as in magnetic resonance force microscopy 38 . 
Reducing the cantilever spring constant by as much as two or- 
ders of magnitude can involve simply increasing the cantilever 
length and clamping it only on a single side. 

Future dynamic AFM measurements will require an actua- 
tion mechanism for driving the cantilever's motion. As an il- 
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FIG. 5: (a) Simulated displacement profile for a disk-double- 
cantilever under application of 50 V to electrodes placed at the sides 
of the cantilevers, (b) Calculated maximum overall displacement and 
displacement along the x-axis as a function of applied voltage. 



lustration of one approach, we consider electrostatic actuation 
through a pair of fixed electrodes that are placed 350 nm to the 
side of the double cantilever geometry (Fig. [5]). Finite element 
modeling shows that stable, steady-state displacements in ex- 
cess of 100 nm can be achieved with an applied voltage under 
50 V. We also note that the displacement is primarily (> 99 %) 
along the x-axis, confirming the effectiveness of the cantilever 
mode engineering described above. In practice, applications 
such as frequency modulation AFM will require much smaller 
displacements, and regardless, the maximum detectable dis- 
placement under the current scheme is w 4 nm. This displace- 
ment level should be achievable for an applied voltage near 
5 V. We can then estimate the performance of this system in 
a frequency modulation AFM scheme using the results from 
Ref. 35, along with the h\ x mode frequency £2m/2tc=2.48 
MHz and linewidth T M /27t=603±4 kHz for the device of Fig. 
Uc). The minimum detectable force F m j n and force gradient 
3F^ n are F min = y/(4kk B TB)/(£l M Q M )=5Ax 10~ 14 N and 
SF^ n = ^{Akk B TB)/(Q.MQMA 2 )=\2x\Q- 5 N/m, where A 
is the cantilever oscillation amplitude (4 nm), and B is the 
measurement bandwidth, taken to be 50 Hz for comparison 
to other experiments 39 . Despite operating in an ambient envi- 
ronment with <2m ~ 4, the estimated F m \ n and &F^ in values are 
competitive with a range of systems operated in an ultra-high- 
vacuum environment. In particular, silicon cantilevers^ with 
k=2 N/m, Q m /2tz=15 kHz, and Q M = l.Ox 10 5 have achieved 
F min = 5.9xl0~ 15 N and 8F^ n = 3.0xl0~ 7 N/m, while 
quartz tuning forks in the qPlus configuration 40 with £=1800 
N/m, Cl M /2%=20 kHz, and Qm = 2.5 xlO 3 have achieved 
F min = 2.1 xl0~ 12 N and 8F^ n = l.lxlO" 2 N/m. More 
recently 39 , ultra- stiff piezoelectric quartz length-extension 
resonators with k=5Ax 10 5 N/m, £1 M /2%=1 MHz, and Q M = 
2.5 xlO 4 have achieved F min = 1.6xl0~ 12 N and 8F^ = 
8.3 xlO -3 N/m. The disk-cantilever system demonstrated 
here operates in an attractive region of parameter space that 
differs from the above sensors, in combining a MHz oscilla- 
tion frequency with a 0.1 N/m to 10 N/m stiffness (with stiff er 
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geometries potentially feasible). 

In summary, we have demonstrated sensitive transduction 
of the motion of a nanoscale cantilever using a high quality 
factor microdisk cavity fabricated on the same device layer. 
Future work will be aimed at understanding the capabilities of 
this system in AFM measurements. This will include mea- 
surements under vacuum to determine ultimate mechanical 
Qm$ of the devices, and to ascertain whether effects such as 
optical cooling and regenerative oscillations^ are accessible. 
Functional devices for AFM will be fabricated to expose the 
probe tips to allow close proximity to other surfaces, and will 
be fully integrated systems combining electrostatic or optical 
actuation with on-chip resonators and waveguides^. 
Methods 

Device Fabrication Devices were created in a silicon-on- 
insulator wafer with a 260 nm thick device layer, 1 jum 
thick buried oxide layer, and specified device layer resis- 
tivity of 13.5-22.5 ohm-cm (p-type). Fabrication steps in- 
cluded electron-beam lithography of a 400 nm-thick positive- 
tone resist, an SF6/C4F8 inductively-coupled plasma reac- 
tive ion etch through the silicon device layer, a stabilized 
H2SO4/H2O2 etch to remove the remnant resist and other or- 
ganic materials, an HF wet etch to undercut the devices and 
release the cantilevers, and a critical point dry to finish the 
processing. The etch time required to go through the sili- 
con device layer is a function of cantilever-disk gap, with an 
« 30% increase in etch time required for G = 50 nm devices 
relative to G = 200 nm devices. 

Device Simulation Mechanical eigenfrequencies and eigen- 
modes of the cantilever and disk were studied using a com- 
mercial finite element software package. Silicon was mod- 
eled as an elastic cubic material using three independent elas- 
tic constants 42 with (100) orientation, and clamped boundary 
conditions were assumed at the cantilever ends. Mesh refine- 
ment studies indicate that numerical errors are below the un- 
certainty resulting from imperfect knowledge of the cantilever 
geometry, which is generally a few percent of the reported 
values. For the reported mode frequencies £Im and effective 
masses m, zero residual stress was assumed. In a separate nu- 
merical study, all cantilever mode frequencies were shown to 
be approximately independent (within a few percent) of the 
residual stress for stress values under ±100 MPa. The mode 
stiffness was calculated as k = mQ^. 

Electrostatic actuation was modeled by iteratively solving 
a coupled three-dimensional static-mechanical problem and 
a three-dimensional electrostatic problem. The former fixes 
the elastic properties and clamped boundary conditions at the 
four double-cantilever ends to be the same. The mechani- 
cally fixed electrodes are 260 nm thick, 500 nm wide, and 
3 jum long, and the gap between them and the cantilever is 
350 nm. The same fixed voltage is applied to both electrodes 
(doped silicon is assumed to be a perfect conductor), while 
the cantilever is assumed to be at the ground potential. Given 
the applied voltages and shape of the deformable cantilever 
and fixed electrodes, the electrostatic force densities on all 
cantilever surfaces are calculated using a boundary element 
method. The calculated force densities were then applied as 
boundary conditions and the mechanical problem was solved 



to find the new deformed beam shape. The electrostatic and 
mechanical solvers were iterated until the solution converged 
to a stable value for each applied voltage. Mesh refinement 
studies were conducted on the electrostatic surface mesh to 
ensure numerical accuracy. The microdisk and substrate were 
assumed to be at ground potential and not included in this 
model for simplicity. This is justified because the cantilever- 
electrode gap is much smaller than the distances between the 
electrodes and either the microdisk or substrate. 

Optical eigenfrequencies and eigenmodes of the disk- 
cantilever system were found numerically using a second 
commercial finite element software package. The silicon 
layer was modeled as having an index of refraction n=3A 
surrounded by air (n =1), and both materials were assumed 
lossless and non-magnetic. The model size was chosen to be 
large enough to fully contain the modes studied, with scatter- 
ing boundary conditions on the outside surfaces. A mesh re- 
finement study was conducted to ensure numerical accuracy. 
gOM for the h\ x mechanical mode as a function of the gap G 
was obtained by linearly translating the cantilever with respect 
to the disk along the x-axis from the initial cantilever-disk gap 
G=100 nm. For each value of G between 30 nm and 300 nm 
the cantilever was further deformed using the calculated h\ x 
mode shape. The modal deformations were and ±d, where 
d varied from 2 nm for G = 30 nm to 10 nm for G > 100 nm. 
For each G and deformation the frequencies and Qs for multi- 
ple optical modes were numerically calculated. For each op- 
tical mode and G the derivative of the frequency with respect 
to modal deformation was obtained using the slope of a linear 
fit. In all cases, the gap changes and cantilever deformations 
were implemented by numerically deforming the same orig- 
inal mesh to obtain the desired cantilever shape and position 
before solving the optical eigenvalue problem. 

Device Characterization Devices were characterized using 
a swept-wavelength external cavity tunable diode laser with 
a time-averaged linewidth < 90 MHz and absolute stepped 
wavelength accuracy of ± 1 pm. The wavelength tuning range 
and linearity are calibrated using an acetylene reference cell, 
so that the uncertainty in optical cavity Qs is dominated by 
fits to the data. Light is coupled into and out of the cavities 
using an optical fiber taper waveguide in a N2 -purged envi- 
ronment at atmospheric pressure and room temperature. Cav- 
ity transmission spectra were recorded using a variable gain 
InGaAs photoreceiver with a typical bandwidth of 775 kHz, 
noise equivalent power (NEP) of 1.25 pW/\/Hz, and gain of 
4.5 x 10 4 VAV. RF spectra were recorded using either a MHz 
(DC) to 125 MHz InGaAs photoreceiver (NEP=2.5 pWA/Hz, 
gain=4x 10 4 VAV) or DC to 1.1 GHz InGaAs avalanche pho- 
todiode (NEP=1.6 pW/yllz, gain=1.4xl0 4 VAV) whose out- 
put was sent into a 9kHz to 3.0 GHz electronic spectrum 
analyzer with resolution bandwidth typically set at 30 kHz. 
RF frequencies and linewidths are determined by Lorentzian 
fits to the data, with uncertainties given by the 95 % confi- 
dence intervals of the fit (uncertainties are not written if they 
are smaller than the number of digits to which the value is 
quoted). Optical frequencies and linewidths are determined by 
a least squares fit to the data using a doublet model that takes 
into account both clockwise and counterclockwise whispering 
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gallery modes and their coupling due to backscatteringaii. 
Phase modulator calibration As a consistency check on the 
calibration of displacement sensitivity 25 , we use an electro- 
optic phase modulator (Fig. Eta)) of known modulation depth 
84> and frequency £l mo d to generate a tone in the RF spec- 
trum, at 44 MHz in Fig. He). This modulation peak is 
equivalent to an effective mechanical oscillation amplitude 
Xmod — 8(|)(co mo j / #om)> and can provide a check on x rms , but is 



limited by the accuracy to which goM is known. For Fig. [2 as- 
suming G = 100 nm and that the optical mode used for trans- 
duction is the TE2,45 mode, goM/27T = 0.61 GHz/nm produces 
a value x « 192 pm that is w 20 % greater than x rm> ? = 160 pm. 
A likely source for the discrepancy is imperfect knowledge of 
the gap; for example, a 10 nm decrease in it would completely 
account for the difference between the two values. 
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Supplementary Information 



1. Optical cavity modes and optomechanical coupling 

Finite-element method (FEM) simulations indicate that the 
p-\ and p=2 modes of TE polarization and p=l modes of 
TM polarization have high Qs (> 10 5 ) for sufficiently thin 
cantilevers. Simulation results for W=65 nm cantilevers are 
shown in Fig. ISlT a). Similar simulations for W=100 nm can- 
tilevers indicate a reduction in Q by as much as a factor of 3, 
though it nevertheless remains above 10 5 . As discussed in the 
text, while most fabricated devices have cavity Qs in the range 
of 5 x 10 4 to 2x 10 5 , a few exhibit Qs as high as w 6x 10 5 (Fig. 
ISTlb)). The optically transduced RF spectra in such devices 
often show a strong amplitude for not only the h\ x modes, but 
also h my and v n modes. This suggests some amount of asym- 
metry in the cantilever structure not found in the majority of 
the devices (such as those studied in the main text). 

Generally, the measured optical Qs decrease with decreas- 
ing gap and increasing cantilever width. Smaller gaps can also 
be problematic because the time required to etch through the 
silicon layer goes up as the gap size is reduced, potentially 
leading to mask erosion and a roughening of the disk side- 
walls. The optomechanical coupling #om, on the other hand, 
increases with decreasing gap and increasing cantilever width. 
The calculated #om for p=l modes with a W=100 nm can- 
tilever is shown in Fig. [Sljc), and can be « 25 % larger than 
the values calculated for W=65 nm in Fig. Old). 



2. Hansen- Couillaud polarization spectroscopy 

For future experiments (including AFM applications) it will 
likely be necessary to lock the probe laser to the cavity. This 
can be done by beating the signal exiting the cavity with a 
strong local oscillator (LO), thereby measuring phase fluc- 
tuations due to cantilever motion and giving access to a dis- 
persive signal needed for locking. A particularly convenient 
approach, Hansch-Couillaud polarization spectroscopy as de- 
scribed in Refs. S2 and S3 and shown schematically in Fig. 
ES2la), sets the polarization so that only part of the input field 
couples to the cavity, with the orthogonal polarization serv- 
ing as the LO. The interference signal is analyzed using a 
X/4 waveplate and polarizing beam splitter, whose outputs are 
measured on a 100 MHz balanced photodetector. The differ- 
ence signal produced by scanning the laser over a cavity reso- 
nance is shown in the inset to Fig. [S2lb). Positioning the laser 
on resonance and measuring the RF fluctuations in this signal 
produces the thermal noise spectrum shown in Fig. [S2lb). 



3. Self-induced optical modulation and free carrier effects 

Two-photon absorption is well-known to play an important 
role in silicon nanophotonics S4 , with the subsequent gener- 
ation of phonons and free carriers giving rise to both optical 
dispersion and loss, and with the associated lifetimes affecting 
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FIG. SI: (a) FEM-calculated optical Qs for TE and TM polarized 
modes of a disk-cantilever with W=65 nm. (b) Thermal noise spec- 
trum of a G=100 nm, W=100 nm disk-cantilever device. To the right 
are two optical modes of the structure; the bottom scan is of the mode 
used in transduction, (c) Predicted optomechanical coupling gQM be- 
tween the h\ x cantilever mode and 1st order radial TE and TM modes 
of the microdisk. Points are from a series of simulations for vary- 
ing G, while solid lines are calculated using the perturbation theory 
method described in Ref . S 1 



the speed of devices intended to exploit these effects. In Ref. 
S5, Johnson et ah observed that under sufficiently strong con- 
tinuous wave input, silicon microdisks of similar dimensions 
to those studied in this work exhibited steady-state oscilla- 
tions in their transmitted power. The authors attributed this to 
competing thermal and free-carrier effects, as the dispersion 
in the refractive index caused by the two effects are oppo- 
site in sign (red-shift for thermal, blue-shift for free carriers), 
and as the cavity mode position shifts due to this change in 
refractive index, the circulating power in the cavity changes, 
thereby changing the rate at which heat and free carriers are 
created. Looking in the frequency domain, the RF spectrum 
of the transmitted signal displayed a number of sharp peaks 
with a spacing of a few hundred kHz. We have observed sim- 
ilar phenomena in our bare (no cantilever) microdisks. Fig. 
E2ta) shows both a broad (up to 200 MHz) and zoomed-in 
(up to 15 MHz) spectrum of the transmitted signal from a mi- 
crodisk with a Q w 3x 10 5 mode coupled to by a fiber taper 
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FIG. S2: (a) Schematic of the setup used for Hansch-Couillaud ho- 
modyne spectroscopy. The input polarization into the cavity is set so 
that a small fraction of the signal is coupled into the mode of interest, 
while the remainder acts as a local oscillator, (b) RF spectrum mea- 
sured using Hansch-Couillaud homodyne spectroscopy. Inset shows 
the difference signal produced when the laser is scanned over the 
cavity resonance. 



waveguide with P{ n ^440 juW at 1533.6 nm. A comb of sharp 
peaks is observed in the RF spectrum, with a nearest-neighbor 
spacing that is typically w 3.23 MHz. 

The devices shown in the main body of the text have some- 
what lower optical Qs than the above device, which likely ex- 
plains why similarly sharp RF peaks are not observed at sim- 
ilar input powers. Instead, the RF spectra (Fig. [2a)) look to 
be a superposition of the spectrum due to mechanical oscil- 
lations of the cantilever and/or disk and the spectrum due to 
competing free carrier and thermal effects within the disk, al- 
beit below the threshold for oscillation. To consider this point 
further, we look at the disk-cantilever device of Fig. ISlT b), 
coupling to a TF^ mode with Q w 1.3 x 10 5 (bottom inset of 
Fig. IS2l b)). If we initially restrict P{ n w 60 and vary the 
laser-cavity detuning AX, we generate fig. IS3f b). We see that 
for initial large detunings (AX = — 35 pm), the RF spectrum is 
dominated by the mechanical modes of the cantilever, but as 
the detuning decreases, a background with broad resonances 
is superimposed (AX = — 31 pm) and dominates (AX = — 23 
pm), before eventually the mechanical modes re-appear for 
small enough detunings (AX = — 15 pm). It is believed that the 
broad background is due to the same thermal/free-carrier ef- 
fects seen in Ref. S5 and in the bare disk of Fig. IS3[ a). Indeed, 
if P[ n is increased to a few hundred juW, a qualitatively similar 
RF spectrum (Fig. IS3[ c)) is observed - here the mechanical 
modes are completely dominated by the thermal/free-carrier 
effects. 

Quantitative modeling of this behavior can be accomplished 
in a manner similar to that of Ref. S6, where the equations of 
bare optomechanics (evolution of the intracavity optical field 
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FIG. S3: (a) Broad and zoomed-in RF spectrum of a bare microdisk 
(no cantilever) with Pi n = 440 juW coupled to a Q « 3xl0 5 cav- 
ity mode, (b) RF spectra from the cantilever-microdisk system of 
Fig. [Sl]as a function of laser-cavity detuning AX with Pi n = 60 jjW 
coupled to a Q « 1.3 xlO 5 cavity mode, (c) RF spectrum of the 
cantilever-microdisk system with Pi n « 1400 juW. Inset is a zoomed 
in high-resolution scan of a portion of the spectrum. 



amplitude and mechanical position) were augmented by an 
equation for the cavity temperature increase. We now have to 
add a fourth differential equation, to account for the change 
in free carrier population. Following the treatment of thermal 
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and free-carrier terms presented in Ref. S5, we have: 
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where for simplicity we have assumed a single-mode cav- 
ity - a more detailed treatment would include both modes 
of the microdisk and their coupling via backscattering. The 
first equation describes the intracavity field amplitude a(t), 
where the first term on the right is its decay due to intrin- 
sic and waveguide loss T, two-photon absorption ((Xtpa) and 
free-carrier absorption (PfcaX while the second term includes 
the laser-cavity detuning 8co c and dispersion due to optome- 
chanical coupling (£om), thermo-optic tuning (g t h), and free- 
carrier dispersion (g/ c ). The second equation describes the 
mechanical motion x(t) with frequency £Im and damping Tm 
and driven by the coupling to the optical field. The third equa- 
tion describes the cavity temperature change AT(t), where the 
cavity has a heat capacity c t h, the temperature decays with a 
rate y^, and is generated by linear absorption (T^ is the por- 
tion of total optical loss that contributes), two-photon absorp- 
tion (octpa), and free-carrier absorption (Pfca). Finally, the 
fourth equation describes the modal free carrier population 
N(t), which decays at a rate y^ c and is generated in propor- 
tion to the square of intracavity energy, with proportionality 
Xfca- The various coefficients in the above equations are de- 
scribed in detail in Ref. S5, and are a combination of physical 
properties such as the two-photon absorption coefficient of sil- 



icon and the absorption cross-section for free carriers, as well 
as cavity mode properties such as its group index and differ- 
ent confinement factors and modal volumes weighted accord- 
ing to the electric-field dependence of the given process (e.g., 
two-photon absorption or free-carrier absorption). 

An analysis of the above equations will produce correc- 
tion terms to the mechanical frequency £Im and linewidth Fm, 
and may help provide a better understanding of the power- 
dependent RF spectra shown in the main text (Fig. [3]). For ex- 
ample, Eichenfield et.al S6 determined that in their SiN x pho- 
tonic crystal nanobeam devices, heating significantly affected 
the linewidth as a function of detuning, but not the mechanical 
frequency. In comparison, in addition to linewidth modifica- 
tion (observation of damping for blue-detuned excitation), the 
frequency dependence on detuning for our devices (Fig. Oc)- 
(d)) does appear to show some effect, in that the shape of the 
curves near zero-detuning is not nearly as sharply- sloped as 
the equations of bare optomechanics predict. 
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